JOURNAL  OF  APPLIED  PHYSICS  105,  023710  (2009) 


Effects  of  electron  scattering  at  metal-nonmetal  interfaces 
on  electron-phonon  equilibration  in  gold  films 

Patrick  E.  Hopkins, a)  Jared  L.  Kassebaum,  and  Pamela  M.  Norris 

Department  of  Mechanical  and  Aerospace  Engineering,  University  of  Virginia,  P.O.  Box  400746, 
Charlottesville,  Virginia  22904-4746,  USA 

(Received  8  October  2008;  accepted  8  December  2008;  published  online  26  January  2009) 

Electron  scattering  at  interfaces  between  metals  and  dielectrics  is  a  major  concern  in  thermal 
boundary  conductance  studies.  This  aspect  of  energy  transfer  has  been  extensively  studied  and 
modeled  on  long  time  scales  when  the  electrons  and  phonons  are  in  equilibrium  in  the  metal  him. 
However,  there  are  conflicting  results  concerning  electron-interface  scattering  and  energy  transfer  in 
the  event  of  an  electron-phonon  nonequilibrium,  specifically,  how  this  mode  of  energy  transfer 
affects  the  electron  cooling  during  electron-phonon  nonequilibration.  Transient  thermoreflectance 
(TTR)  experiments  utilizing  ultrashort  pulsed  laser  systems  can  resolve  this  electron-phonon 
nonequilibrium,  and  the  thermophysical  property  relating  rate  of  equilibration  to  electron-phonon 
scattering  events  G  can  be  quantified.  In  this  work,  G  in  Au  films  of  varying  thicknesses  are 
measured  with  the  TTR  technique.  At  large  fluences  (which  result  in  high  electron  temperatures),  the 
measured  G  is  much  larger  than  predicted  from  traditional  models.  This  increase  in  G  increases  as 
the  him  thickness  decreases  and  shows  a  substrate  dependency,  with  larger  values  of  G  measured  on 
more  conductive  substrates.  The  data  suggest  that  in  a  highly  nonequilibrium  system,  there  could  be 
some  thermal  energy  lost  to  the  underlying  substrate,  which  can  affect  G.  ©  2009  American 
Institute  of  Physics.  [DOI:  10.1063/1.3068476] 


I.  INTRODUCTION 

The  occurrence  of  electron-phonon  nonequilibrium  in 
metal  films  is  an  important  consideration  in  many  nanoappli- 
cations.  For  example,  a  greater  understanding  of  electron- 
phonon  scattering  and  subsequent  energy  transport  has  made 
it  possible  for  microelectronic  engineers  to  develop  held  ef¬ 
fect  transistors  that  achieve  high  radio-frequency  (rf)  power 
levels  in  microwave  radar  and  communications  transmitter 
applications.1  However,  with  these  high  rf  power  levels  come 
extremely  large  thermal  huxes  (>1  kW  cm-2)  that  can  in¬ 
hibit  sufficient  power  dissipation  away  from  active/heat  gen¬ 
eration  regions  causing  self-heating,  increased  operation 
temperatures,  and  thermal  cycling  that  degrade  device  gain 
and  efficiency."  With  the  continued  size  reduction  and  in¬ 
creased  operation  frequencies  envisioned  in  these  devices, 
heat  dissipation  is  becoming  a  growing  challenge  due  to 
large  resistances  from  electron-phonon  coupling.1'3'4  An 
added  challenge  is  dealing  with  the  thermal  resistance  from 
the  electron-phonon  nonequilibrium  at  the  active  layer/ 
substrate  interface.  This  electron-phonon  nonequilibrium  is 
also  critical  in,  for  example,  advancement  of  ablation  and 
laser  machining  of  materials,  7  understanding  spin  dynamics 
in  magnetic  materials,  ^  and  further  development  of  ul¬ 
trashort  pulsed  laser  systems.1 '' 14 

The  thermal  relaxation  between  the  electron  and  phonon 
systems  is  efficiently  observed  with  ultrashort  pulsed  laser 
techniques.15  Electron-phonon  nonequilibrium  resulting  from 
pulsed  laser  heating  can  be  divided  into  three  characteristic 
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time  intervals.16'17  Consider  an  ultrashort  laser  pulse  that  is 
incident  on  the  surface  of  a  solid.  The  earliest  of  the  time 
intervals,  the  length  of  which  is  termed  as  the  relaxation  time 
ree,  is  typically  of  10-100  fs  for  metals.  This  time  repre¬ 
sents  the  time  it  takes  for  the  excited  electrons  to  relax  into  a 
Fermi  distribution  through  electron-electron  (e-e)  collisions. 
These  e-e  collisions  dominate  electron-phonon  ( e-p )  colli¬ 
sions  during  this  time  interval.  Ballistic  transport  of  the  elec¬ 
trons  also  occurs  during  this  time  and  the  depth  to  which  the 
electrons  ballistically  travel  is  significantly  larger  than  the 
optical  penetration  depth  in  s-  and  /7-band  metals.16'16  Once 
equilibrium  is  achieved  within  the  electron  system,  the 
higher  temperature  electrons  transmit  energy  to  the  lattice 
through  e-p  scattering  processes  as  the  electrons  conduct  en¬ 
ergy  deeper  into  the  film  away  from  the  thermally  excited 
region."0'"1  The  e-p  interactions  eventually  lead  to  the  two 
subsystems  reaching  an  equilibrium  temperature  within  a 
time  determined  by  the  specific  heats  of  the  systems  and  the 
electron-phonon  coupling  factor.13  This  electron-phonon  re¬ 
laxation  time  Tep  is  typically  on  the  order  of  1  ps  for  metals 
and  is  inversely  related  to  the  electron-phonon  coupling  fac- 
tor  G,  which  is  typically  on  the  order  of 
1016— 1017  W  m-3  K_1  for  metals.16  Once  e-p  equilibrium  is 
achieved,  thermal  transport  is  accurately  described  by  the 
Fourier  law  as  the  thermal  energy  is  transmitted  deeper  into 
the  film  at  a  rate  proportional  to  the  thermal  conductivity  of 
the  material. 

Electron-phonon  relaxation  at  metal-nonmetal  surfaces 
has  been  examined  by  several  groups,  mainly  in  imbedded- 
metal-nanoparticle  geometries,  with  conflicting  results.  For 
example,  Arbouet  et  al.  found  that  electron-phonon  relax¬ 
ation  time  decreased  (i.e.,  electron-phonon  power  transfer, 
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described  by  G,  increased)  as  nanoparticle  diameter  de- 
creased.  Hodak  et  al.  extensively  measured  electron- 
phonon  relaxation  in  Au  nanoparticles  and  determined  that 
the  diameter  of  the  nanoparticle  had  no  effect  on  overall 
electron-phonon  coupling.  Hodak’s  analysis  showed,  how¬ 
ever,  that  electron-acoustic  surface  mode  interactions  in¬ 
creased  with  a  decrease  in  particle  diameter,  but  the  increase 
was  minimal  compared  to  the  measured  electron-phonon 
coupling  in  the  particle. 

There  have  been  only  a  small  number  of  studies  looking 
at  electron-phonon  relaxation  around  interfaces  in  thin  film 
systems.  These  studies  avoid  complications  due  to  nanopar¬ 
ticle  geometries  (i.e.,  capillary  modes)  on  determining  the 
electron-phonon-interfacial  interactions,  and  still,  groups  re¬ 
port  conflicting  results.  Hohlfeld  et  al.16  studied  G  in  Au  of 
various  thicknesses  for  energies  at  and  near  the  interband 
transition  threshold  (ITT)  at  two  different  incident  laser  flu- 
ences,  —2.2  and  29  J  m-2,  and  measured  the  same  G  values 
as  those  measured  during  low-fluence  intraband  excitations 
by  Smith  and  Norris,26  -2.2  X  1016  W  m"3  K-1.  Hohlfeld’s 
study  determined  that  at  the  specific  probe  energies  around 
the  ITT,  no  electron  temperature  or  film  thickness  depen¬ 
dence  was  observed  in  G  measurements  in  Au  films.  In  ad¬ 
dition,  this  study  verified  an  enhanced  thermal  penetration 
depth  in  Au  due  to  the  large  ballistic  penetration  depth, 
which  is  related  to  the  large  electron  mean  free  path  charac¬ 
teristic  of  the  noble  metals.16'19  This  observation  allows  for 
simplifying  assumptions  in  the  determination  of  G  since 
large  electron  ballistic  penetration  depths  presumably  allow 
for  G  to  be  studied  independent  of  any  other  variables  in 
relatively  thick  films  due  to  the  lack  of  competing  thermal 
diffusion  processes  which  would  complicate  data  analysis. 
Although  Hohlfeld  et  al. 16  did  not  show  a  thickness  or  elec¬ 
tron  temperature  dependence  on  electron  phonon  coupling 
factor  in  Au  films,  Bosco  et  al.  showed  that  e-p  thermali- 
zation  time  can  be  affected  by  the  substrate  thermal  conduc¬ 
tivity.  Bosco’s  study  examined  electron-phonon  relaxation 
times  in  ultrathin  NiFe  films  on  NiO  and  Si  substrates  at  a 
wide  range  of  incident  laser  fluences  and  found  that  rep  in¬ 
creased  with  an  increase  in  electron  temperature.  In  addition, 
this  study  showed  a  substrate  dependence  in  rep  which  can 
be  related  to  the  effusivity  of  the  substrate.  Hopkins  and 
Norris  showed  an  opposite  temperature  dependence  in  reB 
than  Bosco  et  al.  through  electron-phonon  coupling  factor 
measurements  in  20  nm  Au  Alms.  This  substrate  dependence 
was  not  observed  in  the  data  of  Hohlfeld  et  al. 16  most  likely 
because  the  probe  energy  used  during  measurements  of 
Hohlfeld  et  al.16  was  near  the  ITT  energy.  The  thermoreAec- 
tance  signal  around  the  ITT  is  highly  nonlinear  and  changes 
in  optical  reAectance  associated  with  electron  transitions 
around  the  ITT  have  been  shown  to  inAuence  electron- 
phonon  thermalization  measurements.  The  different  trends 
observed  by  Hopkins  and  Norris"  and  Bosco  et  al.  could 
be  due  to  the  different  electronic  band  structures  of  Au  and 
NiFe  which  lead  to  opposite  electron-phonon  coupling  factor 
temperature  dependencies  due  the  different  energies  inducing 

90 

i/-band  to  sip- band  interband  transitions."  Care  must  be  ex¬ 
ercised  to  take  into  account  the  effects  of  these  transitions 
and  the  resulting  reAectivity.  This  study  examines  electron- 


phonon  coupling  factors  in  Au  Alms  in  a  speciAc  transition 
regime  to  conAdently  extract  the  physics  driving  electron- 
interface  and  electron-phonon  scattering  during  electron- 
phonon  equilibration  after  short  pulsed  laser  heating. 

This  paper  examines  the  electron-phonon  coupling  factor 
of  Au  Alms  of  various  thicknesses  on  Si  and  glass  substrates 
measured  with  the  pump-probe  transient  thermoreAectance 
(TTR)  technique  with  a  range  of  laser  Auences.  The  pump- 
probe  measurements  are  taken  at  800  nm  to  only  induce 
intraband  transitions  in  the  Au  Alms,  the  temperature  depen¬ 
dency  of  the  reAectivity  of  which  is  well  known."6  In  addi¬ 
tion,  the  range  of  laser  Auences  is  chosen  to  keep  the  electron 
temperature  in  a  regime  that  does  not  cause  any  i/-band  ther- 
mal  excitations.  A  substrate  dependency  is  observed  in  the 
measurements,  which  is  attributed  to  electron-interface  scat¬ 
tering  during  electron-phonon  thermalization.  This  electron 
interface  scattering,  which  drives  thermal  boundary  conduc¬ 
tance  from  the  high  temperature  electron  system,  is  quanti- 
Aed  by  comparing  the  measured  electron-phonon  coupling 
factor  to  that  predicted  from  traditional  models.  This  electron 
thermal  boundary  conductance  during  electron-phonon  non¬ 
equilibrium  is  directly  measured  by  Atting  TTR  data  to  the 
three  temperature  model. 

II.  THERMAL  AND  THERMOREFLECTANCE  MODELS 
A.  Thermal  model 

The  nonequilibrium  temperature  induced  in  the  TTR  ex¬ 
periments  can  be  accurately  predicted  with  the  two  tempera¬ 
ture  model  (TTM). 18,30  This  model,  given  by  Eqs.  (1)  and 
(2),  describes  the  rate  of  energy  exchange  between  the  elec¬ 
trons  and  phonons  in  a  metallic  Aim 

dT  d  l  dT  \ 

Ce(Te)-f  =  -  urjp)7  -  G(Te  -  Tp)  +  S(x,t), 
dt  dx\  dx  / 

(1) 

Cpd^  =  G{Te-Tp),  (2) 

where  Ce  is  the  temperature  dependent  electronic  heat  capac¬ 
ity  deAned  as  a  product  of  the  Sommerfeld  constant  y  and 
the  electron  temperature  Te,il  G  is  the  electron-phonon  cou- 
pling  factor,  "  Cp  is  the  lattice  heat  capacity  taken  at  300  K. 
which  is  assumed  to  be  temperature  independent  due  to  mini¬ 
mal  lattice  heating  during  the  few  picoseconds  after  laser 
heating,  and  ke{f(Te,Tp)  is  the  effective  electron  thermal  con¬ 
ductivity  that  accounts  for  e-e  and  e-p  scatterings,  deAned  as 
kcff( Te  ,T{)=  keq[Te /  Tp\ . U  This  expression  for  conductivity  is 
valid  for  relatively  low  electron  temperatures  and  is  com¬ 
monly  used  in  TTM  calculations.6’15’26,28’32'35’36  However, 
when  electron  temperatures  are  driven  above  ~500  K,  yet 
remain  much  lower  than  the  Fermi  temperature  so  the  elec¬ 
tron  gas  can  still  be  considered  degenerate,  electron-electron 
scattering  events  will  increase,  which  will  decrease  the  con¬ 
ductivity  from  what  is  predicted  from  the  previous  simple 
expression  for  kefi.  Therefore,  for  temperatures  much  lower 
than  the  Fermi  temperature,  the  more  general  expression  for 
keff  derived  from  kinetic  theory  should  be  used.  This  expres- 
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sion  is  given  by'  ’  Kn  =v2FyTel{3(AeeT;+BepTp)},  where  vF  is 
the  Fermi  velocity  and  Aee  and  Be„  are  electron-electron  and 

5  31  37  38 

electron-phonon  scattering  constants,  respectively."  '  ’ 
The  energy  transport  processes  modeled  by  the  TTM  are 
assumed  one  dimensional  when  the  heated  area  is  much 
greater  than  the  film  thickness.  The  laser  source  term  S(x,t) 
describes  the  radiation  energy  absorbed  by  the  electron  sys¬ 
tem  and  is  given  by15 


S(x,t)=  0.94- 


(1  -R) 


*  pulse*1 


.R\  1  -  exp 


XJ  exp 


--2.77 

R 


d 
R 

t-  2 1, 


\2 

pulse  I 
^pulse  I  . 


(3) 


where  the  film  surface  reflectivity  R,  film  thickness  d,  and 
the  energy  penetration  depth  P  are  material  properties,  while 
the  fluence  J  and  the  pulse  duration  rpulse  are  parameters  of 
the  incident  laser  pulse.  Here,  P  can  be  used  to  describe  the 
sum  of  the  radiation  and  ballistic  penetration  depths.  The 
(1  -exp[(//.P])  expression  ensures  that  all  the  laser  energy 
absorbed  by  the  electron  system  is  contained  in  the  electron 
system  before  electron  thermalization,  and  not  lost  to  the 
underlying  substrate  in  the  event  that  d<P.  This  term  essen¬ 
tially  imposes  insulating  boundary  conditions  to  the  electron 
system  when  t<  ree,  a  time  regime  that  the  TTM  is  not  suited 
to  specifically  analyze.  Equations  (l)-(3)  are  solved  using 

the  Crank-Nicolson  method  and  insulative  boundary  condi- 

22 

tions  are  typically  appropriate. 


B.  Thermoreflectance  models 

In  TTR  experiments,  it  is  the  change  in  reflectivity  A R 
resulting  from  a  change  in  temperature  in  the  sample  that  is 
measured.  The  change  in  reflectance  of  a  metal  can  be  re¬ 
lated  to  the  change  in  temperature  through  the  change  in  the 
complex  dielectric  function  As(&>,  AT)  =  Ae^w,  AT) 
+  /Ae2(w, AT),  where  Ae^to.AT)  and  Ae2(w,Ar)  are  the 
changes  in  the  real  and  imaginary  parts  of  the  complex  di¬ 
electric  function,  respectively.39  For  ultrashort  pulses  (ree 
<fpuls e<Tep)’  the  resulting  reflectance,  and  thus  the  dielec¬ 
tric  function,  is  a  function  of  both  A  Te  and 
ATP:  Ae(u) ,  ATe ,  ATP)  =  Ae  i(u) ,  ATe ,  ATP) 

+  iAe2((o,ATe,ATp).  For  small  changes  in  temperature 
(ATe,AT„  <  150  K),  Ae  can  be  expressed  as  a  linear  change 
in  A Te  and  AT^,15’40'41 


(?£l  (tel 

Ae  =  Aei  +  i  As2  =  — ATe  + - ATn 

dTe  dTp  p 


(te2  (te i 

—AT,  +  —AT,. 


dT„ 


dT„ 


(4) 


where  the  functional  dependence  of  e  is  dropped  for  clarity. 
The  change  in  reflectance  in  TTR  experiments  can  be  related 
to  the  dielectric  functions  by4' 


A  R 
R 


l 

R 


dR  dR 

- Aei  + - Ae2 

(te!  ds2 


(5) 


By  plugging  Eq.  (4)  into  Eq.  (5),  the  reflectance  can  be  sim¬ 
plified  to  the  familiar  thermoreflectance  equation13'40 

A  R 

—  =aATe  +  bATp,  (6) 

K 


where  a^dRI  dTe  and  b<*dR/dTp.  This  reflectance  model  di¬ 
rectly  relates  the  change  in  reflectance  in  TTR  experiments 
to  the  electron  and  phonon  temperatures  in  the  TTM.  The 
coefficient  a  is  determined  by  fitting  the  peak  of  the  electron 
temperature  profile  in  the  TTM  to  the  peak  of  the  TTR  data, 
and  the  coefficient  b  is  determined  by  fitting  the  lattice  tem¬ 
perature  profile  at  some  time  after  the  electron  peak  when  the 
electron  and  phonon  systems  are  in  equilibrium.15'40  Once 
the  change  in  temperatures  of  the  electron  and  lattice  systems 
are  extracted  from  the  experimental  TTR  data,  this  informa¬ 
tion  can  be  used  to  determine  thermophysical  properties  such 
as  the  electron-phonon  coupling  factor  by  fitting  the  TTM  to 
the  TTR  data  using  G  as  a  free  parameter  to  minimize 
error.  9  This  is  discussed  in  more  detail  in  Sec.  III. 

The  thermoreflectance  model  presented  in  Eq.  (6)  is  de¬ 
pendent  on  the  fact  that  A Te  and  A Tp<  150  K.  This  is  easy 
to  ensure  in  metals  with  relatively  high  Sommerfeld  con¬ 
stants,  but  in  metals  with  low  Sommerfeld  constants,  such  as 
the  noble  metals,  very  low  incident  laser  fluences  will  lead  to 
large  changes  in  electron  temperature  yielding  Eq.  (6)  inap¬ 
plicable.  However,  noble  metals  have  a  very  distinct,  high 
energy  ITT  since  the  ,v-band/(/-band  crossing  is  significantly 
lower  than  the  Fermi  level.  The  lowest  energy  (/-band  to 
available  Fermi  level  transition  is  very  large  for  Cu  (2.15 
eV),  Au  (2.4  eV),  and  Ag  (4  eV),  making  it  easy  to  isolate 
the  effects  of  the  interband  and  intraband  transitions  on 
thermoreflectance.4'  Transition  metals,  in  contrast,  pose  a 
greater  problem  in  isolating  the  various  transitions  and  deter¬ 
mining  their  dependencies  on  G  through  TTR  since  charac¬ 
teristically,  transition  metal  band  structures  are  shifted  from 
the  noble  metals  so  the  ,v-band/(/-band  crossing  is  at  or  above 
the  Fermi  level.  This  produces  several  allowable  low  energy 
(/-band  ((/rband)  to  s-band  transitions  making  it  difficult  to 
examine  the  relative  effects  of  inter-  and  intraband  transi¬ 
tions.  For  example,  Cr  has  Fermi  transitions  at  0.8  (ITT),  1.0, 
1.4,  and  1.6  eV,  W  has  transitions  at  0.85  (ITT),  1.6,  and  1.75 
eV,44  and  Ni  has  transitions  at  0.25  (ITT),  0.4,  and  1.3  eV.45 

In  terms  of  the  thermoreflectance  model,  this  high  en¬ 
ergy  ITT  in  noble  metals  can  lead  to  significant  simplifica¬ 
tions.  This  is  apparent  in  the  nature  of  the  complex  dielectric 
function  e.  It  has  been  shown  that  the  change  in  the  complex 
dielectric  function  can  be  expressed  as  Ae=Ae^+Ae^ 
which  separates  the  intraband  effects  [free  electron  effects 
denoted  by  the  superscript  (/)]  from  interband  effects 
[bound-electron  effects  denoted  by  the  superscript  {b)\ 
explicitly.46'47  Both  Ae*^  and  As-h>  are  complex  expressions 
with  real  and  imaginary  parts.  Therefore,  As  can  be  rewritten 
in  the  form  of  Eq.  (4)  as  Ae  =  Ae1  +  /Ae2,  where  As,  =  Asi|/1 
+  Ae(1/,)  and  Ae2  =  Ae2  +A sf\ 

Now  consider  an  example  of  a  metal  irradiated  with  pho¬ 
tons  where  the  incident  photon  energy  is  less  than  the  ITT, 
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TABLE  I.  Properties  used  in  TTM  calculations  for  Au  films. 


7 

[W  m-3  K-2] 

k eq 

(W  m"1  K-1) 

G 

(W  m-3  K-1) 

cp 

(J  m-3  K-1) 

^pulse 

(fs) 

nl,Au 

n2M 

nl,Si 

n2,si 

71.4 

317 

2.2  X  1016 

2.5  X  106 

100 

0.18 

5 

3.692 

0.0065 

R 

R 

R 

R 

30  nm 

30  nm 

50  nm 

50  nm 

P„ 

P 

Au/Si 

Au/Si02 

Au/Si 

Au/Si02 

(nm) 

(nm) 

'O.SiOj 

)!2,Si02 

wl,air 

0.84 

0.88 

0.94 

0.95 

12.5 

120 

1.454 

0 

i 

i.e.,  tioj<  ITT.  In  this  case,  the  change  in  reflectance  of  the 
metal  is  dominated  by  the  intraband  transitions  from  the  free 
electrons  (since  there  is  not  enough  energy  in  the  incident 
photons  to  excite  the  bound  electrons  to  undergo  an  inter¬ 
band  transition).  Therefore,  Ae  =  Ae(A  The  change  in  the 
free  electron  dielectric  function  with  temperature  for  this 
metal  is  described  by  the  change  in  the  free  electron,  or 
Drude  model,  given  by48-50 


2  2  2 
(f)  _  j  ^plasma  _  ^  ^plasma  .  03  ''‘plasma 

&>(&>  +  zcoT)  (or  +  co~)  co(co2  +  co2) 


,  (7) 


where  co  is  the  frequency  of  the  incident  photons,  o>plasma  is 
the  plasma  frequency  of  the  free  electrons,  and  coT  is  the 
collisional  frequency  of  free  electrons  excited  by  the  incident 
photons.  In  order  to  relate  the  change  in  the  dielectric  func¬ 
tion,  and  therefore  the  change  in  reflectance,  to  the  change  in 
electron  and  phonon  temperatures,  and  therefore  relate  TTR 
data  to  the  TTM,  the  temperature  dependence  of  coT  must  be 
exploited.  This  collisional  frequency  is  inversely  propor¬ 
tional  to  electron  collisional  time,  and  therefore  Matthies- 
sen’s  rule51  can  be  used  to  express  the  free  electron  colli¬ 
sional  frequency  as38  AcoT=l/AT~AeeAT2+BepATp,  where 
Aee  and  Bep  are  the  material  constants  relating  to  the  tem¬ 
perature  dependencies  of  the  electron-electron  and  electron- 
phonon  collisional  frequencies,  as  previously 
mentioned.5’31,37’38 

Smith  and  Norris-6  exploited  the  relationship  between 
the  metal’s  reflectivity  and  the  dielectric  function  to  derive  a 
thermoreflectance  model  that  explicitly  takes  into  account 
the  change  in  electron  and  phonon  temperatures  during  intra¬ 
band  transitions,  which  is  reviewed  here.  The  reflectivity  of  a 
solid  is  calculated  by  the  well  known  relationship 

(«1  -  l)2  +  «2 
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R  = 


(jl  j  +  1  )-  + 


(8) 


where  nx  and  n2  are  the  real  and  imaginary  parts  of  the  index 
of  refraction  n,  which,  in  the  case  of  intraband  excitations 
only,  is  related  to  the  complex  dielectric  function  by  n=n] 
+  in2='xS'^ •  The  real  and  imaginary  parts  of  n  can  be  found 
in  tables  of  the  optical  properties  of  materials.  -  However, 
given  a  change  in  electron  and  phonon  temperatures,  the 
change  in  n ,  and  n2  can  be  evaluated,  and  the  change  in 
reflectivity  A R  can  be  calculated.  A RIR  is  calculated  by  di¬ 
viding  A R  by  Eq.  (8),  since  Eq.  (8)  gives  the  baseline  reflec¬ 
tivity  at  a  given  temperature,  and  A R  determines  how  much 
this  baseline  reflectance  changes  as  a  result  of  changes  in 
electron  and  phonon  temperatures.  The  fitting  procedure  of 


this  intraband  reflectance  model  is  similar  to  that  of  the  stan¬ 
dard  reflectance  model  in  Eq.  (6).  In  this  case,  however,  the 
change  in  electron  temperature  is  related  to  the  change  in 
reflectance  via  A coT.  The  change  in  phonon  temperature  is 
still  considered  linearly  related  to  the  change  in  reflectance 
as  in  Eq.  (6)  since  the  change  in  lattice  temperature  is  mini¬ 
mal  in  the  cases  considered  here. 

The  intraband  reflectance  model-'1  gives  the  explicit  elec¬ 
tron  and  phonon  temperature  dependence  of  the  change  in 
reflectivity  of  a  film  undergoing  intraband  transitions  valid 
for  all  temperatures.  In  the  event  that  interband  transitions 
are  excited  and  contributing  to  the  thermoreflectance,  Eq.  (6) 
must  be  used,  and  its  small  change  in  temperature  restriction 
must  be  adhered  to.  The  remainder  of  this  work  will  focus  on 
nonequilibrium  heating  in  Au  when  intraband  transitions 
dominate  the  thermoreflectance.  As  previously  mentioned, 
Au  has  been  extensively  studied,  so  the  data  presented  in  this 
work  have  a  large  database  of  previous  studies  for  compari¬ 
son.  In  addition,  the  intraband  reflectance  model  allows  for 
electron-phonon  coupling  in  Au  to  be  studied  at  relatively 
high  fluences  that  induce  changes  in  the  electron  tempera¬ 
tures  of  thousands  of  degrees  Kelvin. 

Although  the  intraband  reflectance  model  is  suited  to 
study  electron-phonon  nonequilibrium  at  high  electron  tem¬ 
peratures,  in  the  event  that  the  film  thickness  is  approxi¬ 
mately  the  same  as  the  optical  penetration  depth,  Pa 
=  \l4im2,  where  A.  is  the  wavelength  of  the  incident  photons, 
the  reflectivity  on  the  surface  of  the  film  is  affected  by  the 
substrate,  and  therefore,  the  standard  reflectivity  equation 
given  by  Eq.  (8)  must  be  modified  to  take  into  account  the 
optical  properties  of  the  substrate.  Consequently  the  intra¬ 
band  reflectance  model  must  be  modified.  The  Au  films  in 
this  study  are  grown  on  Si  and  Si02  substrates,  both  of  which 
have  very  low  values  for  n2,  the  imaginary  part  of  n,  as  given 
in  Table  I.  This  parameter  is  called  the  extinction  coefficient 
and  indicates  the  optical  absorption  of  a  material.  Therefore, 
the  reflectivity  of  the  Au/substrate  samples  examined  in  this 
work  can  be  calculated  with  the  equation  for  reflectivity  of 
an  absorbing  film  on  a  nonabsorbing  substrate,  derived  by 
Abeles.54  Figure  1  shows  the  predicted  A R/R  of  the  thick¬ 
ness  dependent  intraband  reflectance  model  derived  by  Smith 
and  Norris26  for  30  and  50  nm  Au  films  on  Si  and  Si  CL 
substrates  as  a  function  of  electron  temperature  while  assum¬ 
ing  the  phonon  temperature  is  held  constant  at  300  K.  The 
nonlinearity  with  temperature  and  variability  with  film  thick¬ 
ness  are  apparent.  The  inset  of  Fig.  1  shows  the  reflectivity  at 
300  K  for  Au/Si  and  Au/  Si02  as  a  function  of  film  thickness 
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Temperature  [K] 

FIG.  1.  (Color  online)  A R/R  predictions  from  the  thickness  dependent  in¬ 
traband  reflectance  model  derived  by  Smith  and  Norris  (Ref.  26).  The  non¬ 
linearity  and  substrate  dependency  are  apparent  showing  that  at  high  tem¬ 
peratures  and  thin  films,  this  model  must  be  used  as  opposed  to  Eq.  (6)  when 
only  intraband  transitions  are  induced.  The  inset  shows  the  equilibrium  re¬ 
flectance  R  off  the  surface  of  an  Au  film  on  both  Si  and  Si02  substrates  as  a 
function  of  film  thickness  at  a  wavelength  of  800  nm  (1.55  eV).  For  films 
thinner  than  50  nm,  the  reflectance  deviates  by  more  than  10%  from  the  bulk 
reflectance,  indicating  that  for  films  in  this  range  the  absorbed  laser  energy 
in  TTR  experiments  could  be  affected  by  the  underlying  substrate. 

d.  These  calculations  assume  an  incident  wavelength  of  X 
=  800  nm  (1.55  eV),  the  values  for  nx  and  n2  listed  in  Table 
I  for  air,  Au,  Si,  and  Si02,  and  Aee  and  Bep  for  Au  as  1.2 
X107  K~2  s”1  and  1.23  X1011  K_1  s-1,  respectively.37’55’56 

57 

Note  that  in  Au,  Aee  and  Bep  are  temperature  independent.' 
The  formulation  of  the  intraband  thermoreflectance  model26 
using  the  film  thickness  reflectivity  relations  is  suited  to  re¬ 
late  A R  to  any  AT  in  a  film  of  any  thickness  in  the  event  that 
intraband  transitions  dominate  the  TTR  response. 

III.  ELECTRON-PHONON  AND  ELECTRON-INTERFACE 
SCATTERINGS  IN  AU  FILMS 

This  study  presents  the  measurements  of  the  electron- 
phonon  coupling  factor  of  Au  films  of  varying  thicknesses  on 
Si  and  Si02  substrates  with  the  goal  of  examining  potential 
film  thickness  or  fluence  dependency.  Assuming  an  insulated 
film  system,  the  change  in  film  thickness  and  fluence  will 
alter  the  electron  and  phonon  temperatures  in  such  a  way  that 
can  be  predicted  by  theory.  To  test  this,  Au  films  of  30,  40, 
and  50  nm  of  thickness  were  sputter  deposited  on  Si  and 
Si02  substrates  and  G  was  determined  using  the  same  short- 
pulsed  TTR  technique  and  parameters  described  by  Hopkins 
and  Norris.  The  30,  40,  and  50  nm  Au  samples  were  sub¬ 
jected  to  10  and  15.0  J  m-2  incident  laser  fluence.  Seven  to 
ten  TTR  measurement  scans  were  taken  on  each  Au  sample. 

Phase  corrected  TTR  data58  for  the  30  and  50  nm  Au/ 
glass  film  taken  at  15  J  m-2  incident  laser  fluence  are  shown 
in  Fig.  2.  These  data  are  normalized  at  the  peak  reflectance  to 
observe  the  difference  between  the  electron  cooling  profiles 
at  these  different  fluences,  and  /  =  ()  is  defined  as  the  time  at 
maximum  reflectance.  It  is  the  shape  of  the  cooling  profile 
after  the  maximum  electron  temperature  that  is  related  to  the 
rate  of  electron-phonon  equilibration;  therefore,  the  normal¬ 


Pump-probe  delay  time  [ps] 

FIG.  2.  (Color  online)  Phase  corrected  TTR  data  on  30  and  50  nm  Au/glass 
samples  taking  with  15  J  m-2  incident  laser  fluence.  The  data  are  normal¬ 
ized  at  the  peak  reflectance  to  show  the  difference  in  cooling  profiles  of  the 
electron  systems  after  laser  heating.  The  initial  time  (f0)  is  defined  at  the 
time  of  peak  reflectance.  The  inset  shows  the  calculations  of  the  TTM, 
assuming  homogeneous  heating  [Eqs.  (9) — (1 1 )]  fit  to  the  first  2.0  ps  of  the 
data  after  laser  heating.  The  best  fit  resulted  in  different  G  for  the  two  laser 
fluences. 

ization  of  the  data  is  for  clear  comparison  between  the  two 
data  sets.  The  majority  of  the  electron  energy  loss  to  the 
phonons  occurs  within  1-2  ps  after  laser  heating.  Differences 
in  the  data  are  seen  within  2  ps  after  the  maximum  reflec¬ 
tance,  indicating  a  difference  in  the  cooling  rate  of  the  elec¬ 
tron  system  during  electron-phonon  equilibration. 

In  order  to  determine  G  from  the  TTR  measurements, 
A  R/R  from  the  TTR  data  is  related  to  the  change  in  tempera¬ 
ture  predicted  from  the  TTM  with  the  thickness  dependent 
intraband  thermoreflectance  model  discussed  in  the  previous 
section  and  the  thermophysical  parameters  in  Table  I.  The 
energy  penetration  depth  in  a  material  can  be  estimated  by 
P=yke„IG,  which  gives  120  nm  in  Au  at  room 

^  91 

temperature.”  This  represents  the  distance  to  which  the  ex¬ 
cited  carriers  (in  this  case,  electrons)  penetrate  into  the  film 
before  relaxing  into  a  thermal  distribution  ( t<  ree).  This  can 
also  be  thought  of  as  the  distance  an  electron  travels  before 
experiencing  a  collision  that  redistributes  the  electron  energy, 
which  statistically  is  represented  by  the  mean  free  path  A. 
From  kinetic  theory,59  A  =  3 keq/(CpvL)~  115  nm  for  Au,  in 
good  agreement  with  P,  where  vL  is  the  longitudinal  phonon 
velocity  (electron  velocity  during  diffusive  transport)  taken 
as  3300  m  s-1.55  This  thermal  penetration  depth  was  also 
observed  experimentally  in  Au  films.16  The  Au  samples  in 
this  study  are  thinner  than  this  energy  penetration  depth  d 
<P,  and  therefore  the  electrons  penetrate  the  entire  thick¬ 
ness  of  the  film  before  colliding  and  relaxing  into  an  equi¬ 
librium  thermal  Fermi  distribution.  Therefore,  the  Au  films 
are  assumed  to  be  homogeneously  heated  after  ree,  meaning 
the  electron  temperature  difference  between  the  front  and 
rear  surfaces  of  the  film  can  be  assumed  negligible  {dTJ  dx 
=  0).16  Therefore,  the  TTM  can  be  re-expressed  as 
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and  an  appropriate  choice  of  the  him  thickness  d  ensures  that 
all  the  incident  laser  energy  is  absorbed  in  the  electron  sys¬ 
tem.  Using  the  him  thickness  d  as  the  characteristic  length  in 
which  the  incident  fluence  is  absorbed  in  the  him  is  appro¬ 
priate  when  d  <  P  and  assumes  that  some  degree  of  ballistic 
transport  will  be  blocked  and  reflection  at  the  hlm/substrate 
boundary  results  in  a  higher  energy  density  within  the  him.16 
This  calculation  assumes  total  reflection  of  the  electron  en¬ 
ergy  back  into  the  depth  of  the  him  with  spatial  homogene¬ 
ity.  Using  the  exact  him  thickness  d  in  Eq.  (11)  is  a  slightly 
more  exact  expression  than  P(1 -exp[r//P])  in  Eq.  (3)  and 
can  be  used  when  d<P.  Note  that  P(1  -exp[r//P])  is  meant 
to  reduce  to  d  when  d<P. 

With  the  intraband  thickness  dependent  reflectance 
model  discussed  in  Sec.  II,  Eqs.  (9)  and  (10)  are  fit  to  the 
first  2.0  ps  of  data  immediately  following  the  transient  peak 
and  G  is  treated  as  a  free  parameter  to  attain  the  minimum 
error.  This  corresponds  to  the  time  when  the  majority  of  the 
e-p  coupling  occurs;  however,  excellent  agreement  between 
the  model  and  the  data  is  observed  for  several  picoseconds 
after  the  peak.  The  best  fit  G  is  obtained  at  the  minimum 
error  between  the  data  and  the  fit,  typically  less  than  1%  in 
this  study. 

The  inset  of  Fig.  2  shows  the  TTM  fit  to  experimental 
data  with  the  best  fit  values  for  G  as  determined  from  TTR 
data  taken  on  both  the  30  and  50  nm  samples.  These  data 
show  noticeably  different  cooling  profiles  with  the  only 
variation  between  the  measurements  being  the  film  thick¬ 
ness.  On  the  50  nm  Au/glass  sample,  G  was  determined  to  be 
3.2  X  1016  W  m~3  K-1,  yet  the  best-fit  value  of  G  measured 
in  the  30  nm  Au/glass  sample  was  almost  a  factor  of  two 
larger,  5.9  X  1016  W  nr  '  K"1. 

The  increased  laser  fluence  causes  an  increase  in  elec¬ 
tron  temperature  which  creates  a  subsequent  increase  in  G, 
as  predicted  by  Chen  et  al.b{)  The  predicted  maximum  elec¬ 
tron  temperature  on  the  30  nm  Au/glass  sample  associated 
with  the  15  J  rrf 2  incident  fluence  in  Fig.  2  (along  with  the 
experimentally  determined  value  of  the  e-p  coupling  factor, 
5.9  X  1016  W  m~3  K-1)  is  —1943  K.  At  this  temperature,  G 
is  predicted  as  2.5  X  1016  W  m-3  K~*  with  the  model  of 
Chen  et  al.60  assuming  a  room  temperature  G  of  2.2 
X  1016  W  m  3  K-1;60  this  predicted  value  is  significantly 
less  than  the  measured  G. 

Figure  3  shows  the  values  of  G  measured  on  20,  30,  40, 
and  50  nm  Au/glass  and  Au/Si  samples,  as  predicted  by  the 
TTM  best  fit  with  10.0  and  15.0  I  m-2  fluences  as  a  function 
of  film  thickness.  The  20  nm  data  were  taken  from  Hopkins 
and  Norris.  As  the  Au  films  increase  in  thickness,  the  inci¬ 
dent  laser  energy  absorbed  by  the  electron  system  is 
stretched  out  among  a  greater  film  depth,  thereby  decreasing 


FIG.  3.  (Color  online)  Best  fit  G  measured  on  eight  different  Au  samples 
using  10  and  15  J  m-2  incident  laser  fluence.  Associated  with  each  data 
point  is  the  corresponding  maximum  electron  temperature  determined  from 
the  TTM  fit.  As  temperature  decreases,  the  measured  G  also  decreases,  yet 
the  trends  of  G  with  temperature  show  different  dependencies  depending  on 
the  substrate,  showing  evidence  of  substrate  interference  in  electron-phonon 
equilibration.  The  20  nm  data  are  from  Hopkins  and  Norris  (Ref.  27).  The 
inset  recasts  the  measured  G  on  the  different  samples  as  a  function  of  pre¬ 
dicted  maximum  electron  temperature  achieved  after  pump  absorption  pre¬ 
dicted  by  the  TTM.  The  different  trends  between  the  measured  data  and  G 
predicted  with  the  model  of  Chen  et  al.  (Ref.  60)  suggest  that  scattering  at 
the  film/substrate  boundary  is  affecting  G  measurements.  The  difference 
between  the  Au/Si  and  Au/glass  measurements  suggests  that  the  properties 
of  the  substrate  also  affect  G  measurements. 


the  energy  density  in  the  electron  system.  This  results  in  a 
decrease  in  the  measured  G,  yet  still,  there  is  a  discrepancy 
between  G  measured  on  the  different  substrates.  This  could 
be  a  result  of  high  temperature  electrons  at  the  film/substrate 
boundary  losing  energy  through  electron-boundary  scatter¬ 
ing,  thereby  offering  another  mechanism  of  energy  loss  from 
the  electron  system.  In  the  thermally  conductive  substrate  Si, 
the  thermal  energy  is  more  readily  transferred  from  the  Au 
electron  system  to  the  substrate  than  in  the  thermally  insula- 
tive  substrate  glass.  The  error  bars  associated  with  each  data 
point  represent  the  standard  deviation  from  the  multiple  TTR 
scans  taken  at  each  fluence  on  each  sample.  Note  that  the 
deviation  from  the  mean  increases  as  the  samples  get  thicker. 
This  is  due  to  the  fact  that  there  is  less  substrate  interference 
in  the  film  reflectance  as  the  film  gets  thicker,  and  therefore 
less  of  the  incident  fluence  is  absorbed  (Fig.  1).  This  causes 
a  smaller  thermoreflectance  signal  from  the  Au  surface  since 
the  Au  electrons  are  not  reaching  as  high  of  a  temperature, 
and  therefore  the  signal  to  noise  ratio  from  the  Au  thermore¬ 
flectance  data  decreases.  This  caused  a  greater  deviation 
from  the  mean  G  determined  by  the  fitting  procedure  out¬ 
lined  above. 

The  inset  of  Fig.  3  recasts  the  measured  G  and  Gmeas  as 
a  function  of  predicted  maximum  electron  temperature.  As 
the  Au  film  gets  thinner  and  the  magnitude  of  the  electron- 
phonon  nonequilibrium  at  the  film/substrate  interface  in¬ 
creases,  the  difference  between  Gmeas  and  G  predicted  by  the 
model  of  Chen  et  al.60  becomes  more  pronounced.  The  dif¬ 
fering  temperature  trends  and  values  of  Gmeas  between  the 
Au/Si  and  Au/glass  samples,  which  also  differ  from  what  is 
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expected  in  Au  at  the  given  electron  temperature  via  the 
model  of  Chen  et  al.,60  suggest  that  electron-interface  scat¬ 
tering  is  dependent  on  properties  of  the  underlying  substrate. 
If  Gmeas  was  only  dependent  on  electron  and  phonon  film 
temperatures,  then  results  for  the  Au/glass  and  Au/Si  samples 
would  not  diverge  at  high  temperatures.  However,  the  differ¬ 
ing  high  temperature  trends  between  the  samples  indicate 
that  some  properties  of  the  boundary  (which  would  differ 
based  on  differing  properties  between  Si  and  Si02)  could  be 
affecting  the  measurements.  The  rate  at  which  the  G  mea¬ 
surements  increase  with  temperature  is  much  greater  in  the 
Au/Si  sample  than  the  Au/glass  sample,  which  could  be  due 
to  the  differing  thermal  properties  of  the  substrate.  This  is 
conceptually  intuitive  since  a  thermally  conductive  substrate 
like  Si  would  be  more  apt  to  transfer  energy  away  from  the 
electron  system  than  an  insulating  substrate  like  glass.  The 
higher  the  electron  temperature,  the  greater  the  rate  of  energy 
transfer  between  the  electron  system  and  substrate  through 
electron-boundary  scattering,  thereby  increasing  the  apparent 
electron-phonon  coupling  factor  as  determined  with  the  TTR 
technique.  As  previously  mentioned,  these  trends  were  not 
seen  by  Hohlfeld  et  al.16  who  examined  Au  films  at  various 
fluences  around  the  ITT.  However,  this  channel  of  electron- 
substrate  energy  transfer  could  in  fact  still  be  present  in  their 
data  but  the  large  optical  response  associated  with  the  inter¬ 
band  transition  could  mask  any  effects  of  the  free  electron 
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gas  loosing  energy  to  the  underlying  substrate. 

In  order  to  quantify  electron-substrate  scattering  from 
TTR  data,  the  TTM  must  be  modified  to  take  into  account 
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electron-substrate  energy  transfer  hes.  Hopkins  and  Norris" 
used  a  three  temperature  model  to  quantify  this  energy  path¬ 
way  in  20  nm  Au  films  given  by 

Ce(Te)^7  =  -  G(Te  -  Tp)  -  - f(Te  -  T0)  +  S(t) .  (12) 

dt  d 

To  predict  hes,  Eq.  (12)  is  fit  to  the  TTR  data  using  hes  as  a 
free  parameter  until  a  best  fit  G  is  achieved  that  matches  G 
predicted  via  the  model  of  Chen  et  al.60  The  resulting  hes 
determined  from  the  Au/Si  and  Au/glass  TTR  data  is  shown 
in  Fig.  4.  The  measured  hes  increases  over  the  range  of  pre¬ 
dicted  electron  temperatures,  with  hes  increasing  more  rap¬ 
idly  with  temperature  on  the  Si  substrate  than  the  glass  sub¬ 
strate,  indicating  a  substrate  dependence  on  electron-phonon 
equilibration.  However,  the  thermal  conductivity  is  over  two 
orders  of  magnitude  greater  than  the  thermal  conductivity  of 
glass.  In  fact,  glass  is  a  very  poor  thermal  conductor.  There¬ 
fore,  it  is  not  expected  that  the  glass  substrate  will  act  as  a 
heat  sink  that  will  cause  a  significant  amount  of  electron 
energy  to  transfer  out  of  the  film  during  nonequilibrium  heat¬ 
ing.  The  increase  in  hes  observed  in  the  Au/glass  samples  is 
most  likely  due  to  electron-interface  scattering  at  the  film/ 
substrate  boundary,  which  could  be  facilitated  by  impurities 
and  other  defects  introduced  during  the  deposition  process. 
Although  a  disordered  region  is  still  expected  at  the  Au/Si 
interface,  additional  energy  is  conducted  away  from  the  high 
temperature  scattered  electrons  due  to  the  high  thermal  con¬ 
ductivity  of  Si.  Therefore,  hes  measured  on  Au/glass  repre¬ 
sents  the  effects  of  the  high  temperature  Au  electrons  loosing 
energy  by  scattering  at  an  interface  of  two  materials,  and  the 


FIG.  4.  (Color  online)  hes  determined  from  the  TTR  data  analyzed  with  the 
three  temperature  model.  The  continued  increase  in  measured  hes  on  the 
Au/glass  samples  demonstrates  the  temperature  dependence  of  electron- 
interface  scattering  and  how  if  affects  hes.  The  difference  in  the  Au/Si  and 
Au/glass  trends  shows  the  effect  that  the  conductive  Si  substrate  has  on  hes. 
In  the  thinner  samples,  in  which  high  electron  temperatures  are  achieved 
during  laser  heating,  scattering  events  around  the  film/substrate  interface 
cause  an  increase  in  electron-phonon  scattering.  The  effects  of  these  scatter¬ 
ing  events  are  more  pronounced  on  Si  than  glass  due  to  the  thermally  con¬ 
ductive  nature  of  Si,  and  therefore  heat  is  more  readily  removed  from  the  Au 
electron  system. 

difference  between  hes  measured  on  Au/Si  and  on  Au/glass 
represents  the  rate  at  which  electron  energy  is  conducted 
across  the  film/Si  interface.  This  analysis  assumes  that  the 
differences  in  Au  structure  around  the  Au  or  Si  and  glass 
interfaces  do  not  result  in  significantly  different  electron 
scattering  rates.  Also,  note  in  Fig.  4  that  as  the  trends  in  the 
Au/Si  and  Au/glass  data  decrease  to  a  predicted  electron  tem¬ 
perature  close  to  room  temperature,  the  measured  hes  ap¬ 
proach  values  of  10-100  MW  m-2  K_1,  which  is  the  same 
order  of  magnitude  predicted  by  theoretical  models  for 
electron-substrate  phonon  scattering  assuming  a  negligible 
electron-phonon  nonequilibrium  in  the  film.61'6" 

IV.  CONCLUSIONS 

Electron  scattering  at  interfaces  between  metals  and  di¬ 
electrics  is  a  major  concern  in  thermal  boundary  conductance 
studies.  This  aspect  of  energy  transfer  has  been  extensively 
studied  and  modeled  on  long  time  scales  when  the  electrons 
and  phonons  are  in  equilibrium  in  the  metal  film  and  when 
phonon-phonon  interface  scattering  becomes  a  significant 
mode  of  boundary  conductance.  However,  the  effects  of 
electron-substrate  energy  transfer  have  never  been  exten¬ 
sively  considered  in  the  event  of  an  electron-phonon  non¬ 
equilibrium,  specifically,  how  this  mode  of  energy  transfer 
affects  the  electron  cooling  during  equilibration.  In  this 
work,  G  in  Au  films  of  varying  thicknesses  was  measured 
with  the  TTR  technique.  G  is  measured  as  much  larger  than 
predicted  from  traditional  models.  This  increase  in  G  in¬ 
creases  as  the  film  decreases  in  thickness  and  shows  a  sub¬ 
strate  dependency,  where  G  is  measured  as  larger  on  a  more 
conductive  substrate.  To  explain  this  phenomenon,  electron- 
substrate  energy  transfer  during  electron-phonon  nonequilib- 
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rium  is  considered.  The  three  temperature  model  is  used  to 
quantify  electron-substrate  boundary  conductance  obtained 
from  fitting  the  three  temperature  model  to  the  TTR  data. 
The  trends  in  the  observed  hes  suggest  that  electron  equili¬ 
bration  can  increase  by  interfacial  scattering,  which  is  sub¬ 
strate  dependent. 
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